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Abstract 

The isoscalar structure functions xF^ and F2 are measured as functions 
of X averaged over all permissible for the range 6 to 28 GeV of incident 
(anti) neutrino energy. With the measured values of xF^, the value of the 
Gross-Llewellyn Smith sum rule is found to be F3 dx = 2. 13 ±0.38 (stat) ± 
0.26 {sy St). The QCD analysis of xFg provides A^=358±59 MeF. The 
obtained value of the strong interaction constant as{Mz) = 0. 120^*^4 is larger 
than most of the deep inelastic scattering results. 
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The data on deep-inelastic i/(l7)-scattering in a wide region of momentum trans- 
fer provide a reliable basis for precise verification of QCD predictions In this 
paper the data on the xF^, and F2 structure functions (SF) are presented for 
the kinematic region of relatively small momentum transfer 0.55 < < 4.0 GeV^. 
The value of the Gross-Llewellyn Smith (GLS) sum rule [@] and as (Mz) are 
evaluated. 

The data samples were obtained from three independent exposures of the IHEP- 
JINR Neutrino Detector to the wide band neutrino and antineutrino beams 
of the Serpukhov U70 accelerator. The exposure to the antineutrino beam (i/^- 
exposure) was performed at the proton beam energy Ep = 70 GeV, whereas the two 
i/^-exposures were carried out one at Ep = 70 GeV and the other at Ep = 67 GeV. 
The experimental set-up and selection criteria for CC events are discussed in 0. We 
restricted the range of the measurements in to W"^ > 1.7 GeV"^ and in -^2^(77) to 
6 < E^ip) < 28GeV. The final number of events and the mean values of Q^, (Q^) 
for the three samples are given in Table 1. 

The SF were measured as functions of x averaged over all permissible for the 
energy range 6 < E^(p^ < 28 GeV. Events were binned in intervals of x, and values 
of xF^ and F2 were calculated in these intervals. 

The number of z/^ interactions, n'^, and interactions, n'^, in a given bin of x 
is a linear combination of the 'average' values {F2} and {xF^} of the respective 
SF in this bin (we assume invariance under the charge conjugation): 

= a^-{F2} - 6^-{xF3} 

nh = <,2-{F2} + bl2-{xF,}. 

The subscripts 1 and 2 correspond to the i/^-exposures at Ep = 70 GeV and Ep = 
67 GeV^ respectively. The quantities a"''^ and b'^'" are integrals ('flux integrals') of 
products of the differential neutrino (antineutrino) flux (p'^^'^^ (E) and known factors 
depending on the scaling variables x, y as foreseen by the standard form of the 
differential cross-section for deep-inelastic (z7^) -scattering off an isoscalar target: 

a" = N / l-y - + -y^)E(p''(E)dxdydE , 

b'' = N y{l-^)E(l)^{E)dxdydE 

71 J 2 

etc. Here is the number of nucleons in the fiducial volume of the detector and 
the parameter R = {F2 — 2xF\)j2xF\ measures the violation of Callan-Gross 
relation [^. 

The number Uyi^-^ of neutrino (antineutrino) interactions in a given x-bin was 
obtained from the measured number of neutrino (antineutrino) events in this bin 
corrected for acceptance, for smearing effects arising from Fermi motion and mea- 
surement uncertainties, for radiative effects (following the prescription given by De 
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Rujula et al. |^) and for target non-isoscalarity (assuming d^ju^ = 0.5 To 
determine the appropriate correcting factors the Monte-Carlo simulation of the ex- 
perimental set-up has been carried out using the program CATAS p. We used 
the Buras and Gaemers (BEBC) parametrization [l^ for quark distributions. The 
charmed quark content of the nucleon was assumed to be zero. The kinematic sup- 
pression of d — > c and s —>■ c transitions was taken into account assuming slow 
rescaling |11] and the following charmed and strange quark masses: rric = 1.25 GeV, 
rus = 0.25 GeV. Fermi motion of nucleons was simulated according to ||12[. The 



details of the Monte-Carlo simulation of the known features of the experimental 
set-up are discussed in and [p!3| . 

The number of interactions in a given bin of x is subject to kinematic constraints 
imposed by the cuts in the muon momentum (P^ > IGeV/c 0), in the neutrino 
energy (6 < Ei^^ < 28 GeV) and in the invariant mass square of the hadronic 
system {W"^ > 1.1 GeV"^). These were taken into account in the calculation of the 
flux integrals by appropriate modification of the volume of integration. 

The measured values of xFs and F2 are presented in Table 2 and in Figure 1. 

The systematic errors presented come from the uncertainties of the correcting 
factors due to the choice of some input quark distributions in the event simulation 
program CATAS. These systematic uncertainties were estimated by repeating the 
calculation of the SF using by turns the Field-Feynman and GRV [|l^ quark 
distributions. Note that the systematic errors in Table 2 do not include the normal- 
ization error of 4% for F2 and 11% for xF^. These normalization errors originate 



from the uncertainties in the and z/^ flux determination |]l6 



With the values of XF3, the GLS sum rule ( the integral of F3 ) has been esti- 
mated. Over the interval 0.02 < x < 0.65 it was calculated by numerical integration 
of the measured values of xF^ weighted by ij x. The contribution from the regions 
< x < 0.02 and 0.65 < x < 1 was evaluated by integrating over these regions the 
parametrization of xF^ with the values of free parameters obtained from the fit to 
the data at 0.02 < a; < 0.65. Finally we obtained 

/ — = 2.13 ±0.38 (stat) ±0.26 (syst). (1) 
Jo X 

The systematic error quoted is the quadrature sum of ±0.24 due to i/^ and z7^ flux 
uncertainties, and ±0.09 due to the choice of some input quark distributions. In 
accordance with Table 1 we suppose that the measured value (|l]) of the GLS sum 
rule corresponds to the averaged value ~ 1.1 GeV"^. 

The experimental data on the xF^ were compared with the QCD prediction for 
Q^-evolution by the Jacobi polynomials method in the next-to-leading order QCD 
approximation |]l^ |18|, |19|, ^ . Making QCD analysis of the xF^ SF, for the flrst 



step we do not discuss the problem of validity of application of perturbative QCD 
predictions for kinematical region of small as well as the nuclear effects, heavy 
quarks threshold effects and higher order QCD corrections. 

In order to take into account the target mass corrections the Nachtmann mo- 
ments im of F3 and F2 could be expanded in powers of M'^^^^ JQ"^, and retaining 
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only terms of the order M^^^^i JQ"^ one could obtain: 



N(N ^ 1) A/f2 

M3(2)(iV, Q') = M3«g^(iV, Q') + )^ M3«g^(iV + 2, Q'). (2) 

Here M^^^{N, Q^) and M^^^{N, Q^) are the Mellin moments of xF^ and F2: 

M3(iV,Q2) = [' dxx^'-^xFsix^Q^), 
Jo 



M2iN,Q' 



dxx^~^F2{x,Q^), N = 2,3,... 

Jo 



(3) 



The g2_evolution of M^^^{N, Q^) and M^^^{N, Q^) is defined |2|, |3[ by QCD 
and is presented here for the nonsinglet case for simplicity: 



(In 



2\1 <^JV 
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as (Q2 



Hn[qIQ'')m^^''{N,QI), N = 2,3,... (4) 



0, 



Here as{Q'^) is the strong interaction constant, 7}^^^*^ are the nonsinglet leading 
order anomalous dimensions, and the factor {Qo^ Q"^) contains all next-to- 
leading order QCD corrections p4| . 

The unknown coefficients M3{N,Qq) in (Q) could be parametrized as the Mellin 
moments of some function: 



M^^^^{N,Ql) = [' dxx^-^Ax\l-xy{l+-fx), AT = 2, 3,... 
Jo 



(5) 



where the constants A, b, c and 7 should be determined /^from the fit to the 
data. Having at hand the moments (@) - (|^) and following the method discussed in 
O, [13, M we can write the xF^ SF in the form: 



^ max 



xF^--{x,Q')=x%l-xf Y: e;?'^(x)Ecf (a,/3)M3«^^(j + 2, Q^) 



n=0 



j=0 



where 9"'^(x) are the Jacobi polynomials and c^{a,(3) are the coefficients of the 
expansion of 6°'^(a;) in powers of x: 



j=0 

The accuracy of the SF approximation better than 10~^ is achieved for N„ 
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in a wide region of the parameters a and /? fig, |I9 



4 



Using nine Mellin moments for SF reconstruction and taking into account target 
mass corrections we have determined five free parameters A, b, c, 7 and the QCD 
parameter (Table 3). 

For the - dependence of the GLS sum rule we can write the following theo- 
retical expression]^ : 



GLS{Q^ 



l-a.(Q^)/vr + 0(aD-Al^ 



(6) 



where as is the coupling constant in the MS scheme. The general structure of 



the high-twist (HT) term is known from |26]. The evaluation of this term was 
carried out in [g^], ((O)) = 0.33 ± 0.16 GeV^ , and more recently in ((O)) = 

0.53 ± 0.04 GeV^, using the same three-point function QCD sum rules technique. 
In order to estimate the uncertainties due to the HT contribution we included a 
fenomenological termQ h{x) = x in the fitting procedure . The first moment 

of the function h{x) gives some contribution to the GLS sum rule in accordance 
with p8|. The results of the fit with ((O)) = 0.53 ± 0.04 GeV^ are presented in 



Table 3. The value of as{Mz) was calculated for both variants of the fit due to the 
so-called 'matching relation' p9[. We present the GLS sum rule values calculated 
through with = 1 and with the parameters from Table 3. 

We repeated our fit taking into account both the statistical and systematic errors 



(from Table 2 added in quadrature. With HT from |2^, the following estimations 



have been obtained: A^ = 359 ± 71 MeV, GLS = 2.66 . 

In the singlet case the moments of valence quarks, sea quarks and gluons were 
parametrized at Qq in the form: 

(7) 

M^^^{N,Ql) = dxx^-^Ag{l-xy^, AT = 2,3,... 

Keeping in mind the small number of experimental points we fix Ag from the 
momentum sum rule M^'^^(2,Q^) +M^'^^(2,Q^) = 1 . Following the results 
[|T^ of the QCD analysis of F2 at the momentum transfer Q"^ = 5 GeV^ we put 
Asea = 0.17, Csea = 15 and Cg = 9. The other parameters in (|^) as well as A were 
determined from the fit of the data in the leading logarithm QCD approximation 
and were found to be A^ = 2.49 ± 0.311, K = 0.19 ± 0.02, = 2.80 ± 0.05, 
A = (517 ± 17) MeV with = 6.7 for 6 experimental points and Qq = 3 GeV^. 
Only statistical errors were taken into account. 



■^See for higher order QCD corrections to the GLS sum rule. 

^ This shape of h(x) is in quahtative agreement with the theoretical prediction in and 
experimental estimations in for the x values from Table 2. 
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Several comments: 



The values (|l]) and the results on the GLS sum rule in Table 3 are consider- 
ably smaller in comparison to the results of previous measurements. (See the 



summary on the GLS sum rule data in [32| and the latest 3-loop result [p3|. 



The parameter Ajjg is found to be about twice as large as the estimations 



in and |3|]. It is in qualitative agreement with the results of the NLO 
analysis ||3^ of the GLS sum rule in the MS scheme: = 317±23{stat) ± 

99{syst)±62{twist) MeV with HT and A^^ = A35 ±20 (stat) ±87 (syst) MeV 
without HT. 

• The illustrative nature of the QCD fit to the data on F2 should be pointed 
out. The matter is the absence of reliable theoretical predictions for HT con- 
tribution to singlet SF. In spite of this, we obtained the momentum fraction 
carried by quarks in the nucleon, M^^^ {2,Q'^) = 0.46, to be in agreement 
with the previous measurements. 

• The strong interaction constant at the point of Z boson mass is found to be 
higher than most of the deep inelastic scattering results |36|, |3^ . 

• The consideration of the HT contribution decreases and appreciably changes 
the parameters of the fit as well as the GLS sum rule value and as{Mz)- For a 
reliable QCD analysis one must calculate not only the GLS sum rule (A^ = 1) 
but also the higher SF moments (A^ = 2,3,...). Using in addition a 3-loop 
QCD analysis one could expect to improve the estimation of as{Mz)- 

In conclusion, let us stress once more that the QCD analysis of SF is sensitive 
to the HT contribution and in the future it should take into account the nuclear 
effects, heavy quark threshold effects and higher order QCD corrections. We hope 
to improve the accuracy of our estimations by processing the additional data on 
deep inelastic scattering obtained with the IHEP-JINR Neutrino Detector in the 
wide band beams of z/^ and 77^. 
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Beam 
Ep (GeV) 
Final statistics 



Table 1. Summary of the exposures. 







^/^ 


70 


70 


67 


741 


2139 


3848 


1.2 




2.3 
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X 




F2 


stat 


syst 


AF2 




stat 


syst 


.052 


.55 


1.169 


.026 


.047 


.023 


.445 


.044 


.062 


.148 


1.4 


1.097 


.026 


.022 


.022 


.583 


.044 


.017 


.248 


2.2 


.894 


.023 


.018 


.019 


.622 


.038 


.019 


.346 


2.9 


.576 


.016 


.017 


.013 


.556 


.027 


.011 


.447 


3.4 


.390 


.014 


.012 


.009 


.336 


.023 


.007 


.563 


4.0 


.182 


.008 


.004 


.004 


.177 


.012 


.005 



Table 2. The isoscalar structure functions F2 and xF^ obtained on the as- 
sumption of R — 0. The difference AF2 between the values of F2 obtained with 
R — .1 and those obtained with R — is also presented. The bin edges are at 
X = .0, .1, .2, .3, 4, .5, .65. 
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((O))=0 


((0)) = 0.53 H 




2.8 


2.05 


A 


9.28 ± 1.73 


0.90 ± 0.67 


b 


1.06 ± 0.11 


0.31 ± 0.18 


c 


3.22 ± 0.31 


3.64 ± 0.21 


7 


-0.90 ± 0.21 


9.53 ± 5.73 


[MeV] 


417 ± 51 


358 ± 59 


GLS sum rule 


1.59 


2.63 


as{Mz) 


0.123l^ 


0.120t^ 



Table 3. The results of the NLO QCD fit to the xFg SF data for f = 4, Ql = 
3GeV^, Nmax = 12, a = 0.7, /? = 3.0 with the corresponding statistical errors. 
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Figure captions. 



Fig.l. The x- dependence of the isoscalar structure functions ^2(0;) and 0:^3 (x). 
The statistical and systematic errors are added in quadrature, excluding the nor- 
malization error of 4% for F2 and 11% for xF^. The curve is fit of the form 
xF^{x) = Ax''{l — x)'^. The best fit values of free parameters A = 5.36±1.25 (stat), 
b — 0.81 ± 0.10 (stat), c — 3.52 ± 0.26 (stat) were obtained using for each x-hin the 
mean x of the bin as the actual x-point corresponding to the value of the structure 
function obtained. 
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